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We consider a Bose-Einstein condensate (BEC) of -/V atoms of repulsive interaction ~ Uo, in an 
elliptical trap, axially pierced by a Gaussian-intensity laser beam, forming an effective (quasi-2D) 
toroidal trap with minimum at radial distance p — p p . The macroscopic angular momentum states 
^i(p, 9) ~ \fN^i{p)e M for integer / spread up to p < p max ~ (NUo) 1 ^ 4 2> p v . The spreading lowers 
rotational energies, so estimated low metastability barriers can support large I < lmax ~ (NUo) 1 ^ 4 , ^ 
10 for typical parameters. The /-dependent density profile j$i(p)| 2 — |<3>o(p)| 2 is a signature of BEC 
rotation. Results are insensitive to off-axis laser displacements po, for po/pmax -C 1. 

PACS: 03.75Fi, 05.30.Jp, 32.80.Pj 



Bose-Einstein condensates (BEC) of atoms in magnetic 
traps have been the focus of intense recent activity [|l] |l2| . 
Coherence of the BEC wavefunction was demonstrated 
by formation of atomic interference fringes, on switch- 
off of weakly-coupled traps j|. Non-destructive tests 
of phase-coherence, through the Josephson-like effects 
have been proposed [|)-^), including novel self-trapping 
phenomena Other phase-coherence signatures in 

trapped BEC would be of much interest. 

In bulk superfluids /superconductors, equilibrium 
quantized vortices jl3],[l4|], provided classic evidence of 
phase coherence. In multiply-connected axial-hole ge- 
ometries (e.g. thick superconducting cylinders |15| or su- 
perfluids in narrow rings), superflows with wavefunction 
<f>i(6) = \<P\e iW , and |$| uniform Jl| - |hg[ ] maintain integer- 
/ metastable quantized rotational states. BEC I ^ states 
in a He II like uniform-density limit have been considered 
with the BEC rigidly restricted to flat-potential 
regions by (narrow-ring) rigid walls. However, one needs 
to go beyond this 'square- well container' limit to capture 
the characteristic ./V-dependent spread of BEC in a poly- 
nomial trap. The Gross-Pitaevskii equation (CPE) JToj ] 
for a non-uniform BEC includes both (repulsive) atomic 
collisions that spread BEC, and trap potentials, that 
confine it. GPE vortex solutions centered in (simply- 
connected) harmonic traps have been found Jl(],[ll|], but 
could be unstable to outward vortex displacements fll^l . 

The investigation of TV-atom non-uniform metastable 
states, of macroscopic angular momentum Nlh, in 
multiply- connected toroidal traps are thus important as 
possible signatures of BEC phase coherence. Central the- 
oretical issues include the existence of such GPE states; 
the mixing of /-states by nonlinear GPE terms; the effects 
of off-center (azimuthally asymmetric) displacements of 
the toroidal hole; the existence of metastability barriers, 
and the decay-channel limiting of accessible /- values. 

In this Letter, we consider N atoms in the BEC, of 
interatom scattering ~ Uq > in an elliptical trap, 
with an axial 'hole' (or strong potential barrier) drilled 



by an intense off-resonant laser beam, forming an ef- 
fective toroidal trap. The T = GPE is solved for 
a particular quasi-2D toroidal trap, with the Thomas- 
Fermi approximation (TEA) outside an effective 
core, and a polynomial solution inside it. The cylindrical- 
coordinate wavefunctions ^i(p,6,z) ~ \/~N$x(p)e M are 
used as bases for variational states. A dimensionless ex- 
pansion parameter 5, that is small for strongly interact- 
ing and dense trapped BEC, appears naturally. We find 
that the I = state <fr;=o(p) of energy E Q ~ S" 1 peaks at 
the toroidal axis, spreading out to a radial distance p = 
Pmax ~ S^ 1 / 2 . The centrifugal force induces a density 
fall-off \MP)\ 2 - l^=o(p)| 2 ~ l 2 5 3 [HS- 2 ) - p 2 ma J P % 
that is a rotation signature of the I ^ BEC states, of 
rotational energy E t — E ~ / 2 (ln<5)<5 ~ I 2 / Pmax- Off- 
axis displacements ~ po of the toroidal 'hole' are unim- 
portant if pa/ p m ax -C I- The self-interaction £/ | , I'| 2 , 
rather than destroying /-states, provides metastability 
barriers Eb against Z — > Z — 1 reductions. Other decay 
channels, by emission of unit- vortices, or of Bogoliubov 
quasi-particles, have comparable thresholds, Eb ~ 0(1). 
Increasing the BEC spread p ma x lowers the / ^ cost, 
Ei - E ~ l 2 /p 2 max , so ~ 0(1) barriers E B > E t - E 
can still maintain large / values, / < lmax ~ <5 -1 / 2 . 
These quasi-2D trap results, as well as a square- well con- 
tainer limit, are recovered, from a general d-dimensional, 
anharmonic-trap scaling analysis. We now outline argu- 
ments for the above conclusions. 

i) GPE and I -states: The macroscopic condensate 
wavefunction ^(p,9,z) = y/~N&(p,B,z) obeys the GPE 
in cylindrical coordinates: 

" n 2 ( d 2 id id 2 d 2 \ 

2m \dp 2 pdp p 2 d6 2 dz 2 J 
V trap (p, 9) + V trap (z) + U Q \^\ 2 ] * = E% (1) 

where Uq = 4irh 2 a/m and aim) is the atomic scatter- 
ing length (mass) and E is the (single-particle) energy. 
Vtrap(p, 0)(Vtrap(z)) is harmonic with curvature uj\\(lu z ). 
However, Vt rap (z) rises sharply for \z\ > L z /2, as in 
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Fig. [|. The interacting BEC cloud is blocked from ex- 
panding beyond \z\ < L z /2, spreading, for uJ^Pmax ^ 
oj 2 L 2 , only in the p-direction with increasing N. The 
wavefunction is quasi-2D, uniform ~ \j\fL~ z in the z- 
direction. An off-resonant gaussian-profile laser barrier 
of high intensity V c , and narrow width 2a, directed along 
the z-axis, yields an axially symmetric 'doughnut' trap 
V t rap{p,0) -> Vtrapip), where 

Vtrap(p) = \muj\p 2 + V c e- p2 ^\ (2) 
The BEC wavefunction is, with I = 0,1, 2, 

Hp) 



<S>i{p,0,z) -> 



(3) 



Pma * ^i$ ; (p)i 2 + / Pc d^i*r(p)i 2 = i. (7) 



With the first integral dominating, we have to 0(5 2 ), 

E[ = E + \l 2 5 : Eq = 2N6 = {4:N) 1/2 . (8) 

Here, A = \n(p max /p c ), p c ~ a[ln(V c / Eq)} 1 / 2 , so the 
results are only very weakly dependent (A ~ lnp c ~ 
ln[ln[t^/o']]) on the toroidal core region and how it is 
modeled. The interaction pushes the BEC far beyond 
the toroidal axis, to p — p max {l) p p , and with 

Pm = Prnaxi 1 = °)> we have 



We scale lengths (energies) in the V c = harmonic 
trap length ry = (TL/mwu) 1 / 2 (energy = huiu/2), so 

p~ = /°/ r H> 5 ' = v/ r \\i E l = Vc = K/^^ii. Then 

Eq. (hi), in dimensionless form, is 
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pdp J 
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n\Hp)\ 2 ] Hp) 



^trap{p) + 



(4) 



with N tx N defining a central dimensionless parameter, 



5 = N 1/2 = (iV /A0 1/2 ; N = L z /Ac 



(5) 



As shown in Fig. g, the toroidal trap potential Vt rap (p) = 
p 2 + V c e-P 2 / 2 ° 2 of Eq. (g) has a minimum around the 
(circular) toroidal axis at p = p p = \i2a\h\(V c /2a 2 )] 1 / 2 : 
V trap (p) « (p 2 /a 2 )(p - p p ) 2 + F p , where F p = p 2 + 2a 2 . 
A non-interacting BEC in such a minimum has a gaus- 
sian width r p = r\\{a/p p ), i.e., the 2D 'volume' occu- 
pied is ~ r 2 . The interacting BEC spreads (as seen 

later) to p max ~ N 1/4 , with a volume p 2 max > r|. 
Thus S of Eq. (||) is essentially the BEC volume ratio, 
S ~ T'y/Pmarr ^ I- Note that (5 ~ 1 for the 'square-well' 
container fl7| , |l8l , pof , where the BEC cannot expand. 

We have solved Eq. (^) numerically. In the regime 
5 <C 1, and for a sufficiently low- width laser-hole, a <C 
<5~ 2 / 13 , we find analytic results can be obtained. Drop- 
ping derivatives ~ 6 <C 1 in Eq. (Q), we have, for p such 
that p n 



the TFA wavefunction 



Hp) 



> P> Pc> Prnim 

^) ' [1 + [£i ~ 1 2 /P 2 ~ V trap (p)] /E 



1/2 



(0) 



where Ei = E\ — E n is the I ^ energy per particle above 
the ground state and Hp ma x) = Hp min ) = °- The 
wavefunction is zero for p > p max ; it is matched for p < 
p c to a polynomial solution H inside an effective 'core' 
Pc ~ Pmin- HP) P eaks at p = p peak ~ p p + {a 2 /pj)l 2 . 
The energies are determined by normalization, 



The local density difference due to centrifugal effects ~ 
l 2 /p 2 in Eq. (|) is 

\Hp)\ 2 - |l> ; =o(p)| 2 = M 2 S 3 [1 - (P 2 m /2\p% (10) 



with crossing $ at p w p m /V2A. 

Returning to the p < p c series solution, H(p) — 
J2i=o a 2iP 2l+l , we find that the coefficients are ao ~ 
a 4 /V c 2 , a 2 - a 4 /y c , with a 4 ~ a 3 / 2 S. The TFA_$/(_p) falls 
off as ~ {p-p min ) 1/2 , with p mm ~ V2a[hi(Vc/Eo)] 1/2 . 
The wavefunctions and derivatives can be matched at a 
Pc = [{I + 4)/(Z + 7/2)]p ml „ > p mi „, and the $f integral 
in Eq. (0) is ~ alpj ~ a 13 ln(y c 5)5 2 < 1. 

Figure || shows |$;(p)| 2 versus p, for Z = and 10, 
for parameters as given later. The analytic approxima- 
tion (solid line) closely matches the (dashed line) nu- 
merical solution of Eq. (|J). The centrifugal reduction 
and shift of the density peak (~ 10% for / = 10), is 
a rotation signature, possibly detectable by phase con- 
trast imaging Q. The super-current density is j ~ 
h/{im2-Kp)^*-§g'^ = N\^i(p)\ 2 v s s, where the azimuthal 
velocity v s g = (hl/2iTmp), and j vanishes at the origin 
as j ~ p 2l+3 . The laser hole 'pins' the azimuthal velocity 
of average angular momentum (^>i\(fi/i)d/d6\^>i) = Nl%, 
suppressing displacement instabilities p"2| . 

ii) Metastability barriers: The nonlinear term in 
the CPE could cause mixing of Z-states and I — > 
I — 1 decays, but instead, induces barriers. A GPE 
trial function for full superposition is *Sf(p,0)/\/N — 
[C;i>;(p, 9) + C;_i$;_i(p, 6)], with variational constants, 
\Ci\ 2 + |C;_i| 2 = 1. With an external angular velocity 
51 = dun/2 in the hamiltonian H — > H — TiVlL, the free 
energy functional F(^>) (whose variation is Eq. (|l|)) is 
given by, 



F 

N 



= {E P (n)-i p , p (i-\c p \ 2 )}\c p 

p=i,i-i 

+ 4|C ; | 2 |C ; _ 1 | 2 / U -i, 



(11) 



where I u> = ^N((\Hp)\ 2 \H(p)\ 2 )), _with «...)) = 
/ °° dpp d ~ x and d = 2 at present, and Ei(Q) = Ei — lCl. 
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The Z-state (C/, C;_i) = (1, 0) is separated from the I— 1- 
state (0, 1) by a barrier of height F B /N = E B = E — 
E t (Q), larger than the splitting, AE t = (2/ - 1)A<5. The 
'full-overlap' barrier E B = E B i is at \Ci\ 2 = |Cz_i| 2 = 

E B1 = 2E [1 - A(Z 2 + 1) 2 )S 2 } - (A^ - H). (12) 

- - - 1/2 

Thus Bsi ~ E ~ AT' , and the Z-state (1,0) is pre- 
served by metastable barriers for £~i — l 2 X/"vN < E B i 

-1/2 

or I < N . These barriers arise because interatomic 
collisions disfavor the transition-state slowing down of a 
minority fraction (< 1/2) of the rotating flow. 

We now consider a superposition of non-uniform Z- 
states, overlapping only at an interface, lowering the I, I — 
1 interspecies scattering. With vi, vi-i, the fractional vol- 
umes occupied, and vq = (yi + v\—\ — 1) > the overlap 
fraction, we take $ p (p) — > $> p (p)g p (p) where the step- 
function gi-i(p) (gi(p)) is unity for vi-ip 2 rn > p 2 > p 2 c (for 
Pm > P 2 > (1 — v i)Pm) an d zero otherwise. Repeating 
the previous argument with these wavefunctions, we find 
that the intra-species scattering is increased by the con- 
finement to a reduced volume VpPm < Pm, so the TFA en- 

l)/v p . The 



ergies are raised, E p (v) ~ I p ^ p (v) ~ 

- 1/2 

inter-species scattering gives Iij~i(v) ~ i/q-ZV jviVi-\. 
The lowest interface-overlap barrier i?B2 obtains for ap- 
proximately equal volume fractions ^; = f;_x ^ 1/2, 
and for an annular overlap volume of radius p m j\[2 and 
thickness ~ £. Thus = y/2£p m ; E B 2 is a fraction 
of the full-overlap barrier, E B2 ~ (\/2^/ p m )E B i. (Gra- 
dient contributions omitted in TFA would be compara- 
ble). Here £ ~ l/[A>|l> ; ( / 5 p )| 2 ] 1 / 2 - l/^ 72 is a 'heal- 
ing length', so 75^2 = 2y2 ~ O(l), maintaining Z-states 
for Z < Z ma2; 2 = (2\/2/A) 1 / 2 ( 5- 1 / 2 ,- A 1/4 . The exter- 
nal angular velocity induces an Z — ► Z — 1 transition only 
above critical values, |f2| > |f2 c /| = |(i?B — AE{)\. This 
^c; ~ E B 2 estimate exceeds an Ct c i ~ lAE^I estimate 
based on the splitting alone. 

Wedge-shaped Z — 1 regions, expanding azimuthally, 
would have ~ (2£,/p m )E B i interface costs, and similar 

- 1/4 

Imax ~ N scaling. 2D phase-slip saddle-point solu- 
tions |l^,^(| would also provide such barriers. 

iii) Off-center toroidal hole: In practice, the toroidal 
'hole' of Fig. (|l|) could be drilled off the elliptical-trap 
axis, or the laser could fluctuate in profile and posi- 



tion. With a displacement p — (po , 6) of the 



el- 



liptical trap, V t ra P (p) — * Vtrap + P - 2p pcos9, modi- 
fying the TFA solution of Eq. @. Repeating the above 
argument, the BEC cloud becomes anisotropic, extend- 

- 1/2 

ing to p m = p cos 9 + E , but global averages of 
the energy /angular momentum are unaffected, to order 
~ p /p m <C 1. Although CiC\_ Y terms (previously zero 
by e lW orthogonality) now enter in Eq. ( |TT| ) , the barriers 
E B i t 2 to Z-decay are similarly insensitive. 

iv) I -decay channels: An Z-state could decay through 
(rectangular) unit-vortex loops, with straight segments 



of vorticity J z = — 1 reducing core vorticity by I — > I — 1, 
that pushes the J z = 1 segment outwards by a repulsive 

energy^ (Z-l)lnp. With |$(p c )| 2 cx [p c -p m ) ~ 

p c ~ a, the nucleation cost near the core is E B3 ~ 
£ A^|l>;(p c )| 2 ~ a. Then £/ = AZ 2 <S < E B3 accesses this 
decay channel only for Z > l maX 3 ~ sfdN 1 ^ j '\/X. (As 
v s ep~ ~ M/2inn is a constant, i^e should not induce a 
loop-expanding current-drive force fl4}| .) 

Another possible decay channel is a successive reduc- 
tion N — > N — 2 of BEC Z-state atoms, producing surface 
pairs of dissipative quasiparticles jll[] of Bogoliubov ex- 
citation energy 2 A, and angular momentum l±q ]l7| , ^8| , 
that slow down to fall into the Z = ground state. A de- 



tailed analysis |16| for toroidal traps is beyond the scope 



of this paper; however, for £i < A, the quasiparticle 

In other con- 



channel is not accessed for Z < VAN 1 ^ 4 . 
texts |2l| ], quasiparticle damping is small. 

v) Numerical estimates: Parameters chosen are u\\ = 
132rad s _ , (InK), r« — 6.3pm, L z — 25pm, 2a = 
12pm, a = 50A, V c = 63nK, and m = 2 x l(T 26 kg 
(~ sodium). Then p p = 17.5pm, Nq sa 10 3 , and for 

ZV = 10 6 , we have S = 1/Vn = 0.03. Physical magni- 
tudes are then E = 32nK, X — 2,£i= 0.031l 2 nK, p m = 
50pm, £ = 0.8pm, v se (p) < v se (p c ) = 0.12Zmms _1 , Q c i = 
30 Hz, and l maX 2 ~ 10. Thermal excitations activated 
over NE B free energy barriers are exponentially dilute, 
at low temperatures. 

Finally, we present a general TV-scaling argument. 
Consider a d-dimensional toroidal trap that, for distances 
f well away from the core, is Vt rap ~ hu>tf a , with f = r/r t 
and rt = (h/muit) 1 ^ 2 - For a = 2 the trap is harmonic; for 
a — » oo, a square- well container results, Vt rap = 0(oo) for 
r < r t (r > r t ). Then defining N = (NU /hco t )rf, $(r) = 
<S>{r)/r d t l2 ,Ei = Ei/hcut, the TFA wavefunction, $ ; = 
{[Ei -_f a -(l 2 /f 2 )]/N) 1 / 2 is almost flat, |§,(f)| ~ |#o| » 
(Eo/N) 1 / 2 , with BEC spreading until f ~ f m , when 
the confining potential rises to the energy, 75; ~ f^. 
With normalization ((|$;(f)| 2 )) ^ l^ol 2 ^™ = 1, we ob- 



tain r„ 



TV 



l/(t»+d) 



ing length is £ 
integral J; 



Then £ t - 
(iV|$ | 2 )- 1/2 



Z 2 /f 
- N 



2 n , and the heal- 

a /2{a+d) 

^((l*!^)! 3 !*!-!^)! 3 )) defines a full- 
overlap metastability barrier E B1 ~ ~ (^/^m)> 
while for a cZ— 1 dimensional Z, Z — 1 interface of thickness 
£iF B 2/N = E B 2 = (£/f m )E B i. The vortex nucleation 



1,1-1 



threshold is E m - £ f^7 2 /V|$ 



,-(<i-2) 



The Z-state 



energies 5/ exceed these barriers only for Z greater than 

linncf 



A j-(2+a)/2(a+d) , 



1 r T d/2(a+d) 

N . For a 



A v(2+a)/4(a+d) , 



d — 2, we recover our results, with 
lmax2 ~ Z ma;C 3 - A 1/4 . For a ->■ oo,f m = 1, i.e. a re- 
stricted r m = r t , and the familiar He II healing length 
£ ~ (NUq)^ 1 / 2 is recovered, with energies Ei ~ A + Z 2 , 
and barriers - iV,E B 2 ~ N 1/2 ,E B3 - 0(1). The 



Z-bounds are l maxl = N 1/2 ,l max2 ~ TV 1/4 ,Z maa; 3 - 0(1): 
narrow-ring containers |l7],[l8|,^0| would suppress vortex- 



3 



loops and make quasi- ID phase-slips [g0| the lowest decay 
thresholds, Fbi ~ N 3 ^ 2 Uq^ 2 . Barriers vanish in the ideal [l^ 
gas limit, Uq — ► 0. 

Experimental preparation of I ^ states could be 
through stirring normal states by a laser 'paddle-wheel' 
[ l| and cooling below transition; or by rotating the trap 
[22] with Q > fl c i, below T c . Experimental detection 
might be through the / ^ density profiles, or the Sagnac 
effect |Q in large azimuthal-velocity regions around the 
core. Another interesting setup is a double-toroidal trap, 
with Ni = N,h = 1(JV 2 = N,l 2 = 0) for toroid 1(2), 
where the toroids initially are separated by a block- 
ing laser sheet, whose removal could result in inter- trap 
'macroscopic quantum coherence' p3| . 

In conclusion, I ^ macroscopic angular momentum 
states of condensates, with /-dependent density profiles, 
can exist in toroidal traps, maintained by atomic colli- 
sions, and robust to toroidal hole displacements, and to 
estimated decays, for I not too large. The observation of 
these states would be a signature of macroscopic phase 
coherence of trapped Bosc-Einstein condensates. 
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max i 
FIG. 1. Schematic plot of quasi-2D 'doughnut' toroid ge- 
ometry with laser along the z axis. The trap potential rises 
sharply (dashed boundary) at L z . Here p c ,p P , and p ma x 
are the effective core radius, toroidal axis, and extent of the 
shaded BEC cloud respectively. 




FIG. 2. Trap potential (right vertical axis, dash-dot line) 
Vtrap(p) versus cylindrical-coordinate radius, p,(with bars de- 
noting scaling in harmonic trap energy /length). p is the ef- 
fective toroidal trap minimum. Scaled BEC density |4 > !(p)| 2 
(left vertical axis) versus p from analytic (solid line) and nu- 
merical (dashed line) results. Parameters are as in text. 
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